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Abstract 

A linear point design methodology for minimizing the error 
in on-line Kalman filter-based aircraft engine performance 
estimation applications is presented. This technique 

specifically addresses the underdetermined estimation 
problem, where there are more unknown parameters than 
available sensor measurements. A systematic approach is 
applied to produce a model tuning parameter vector of 
appropriate dimension to enable estimation by a Kalman filter, 
while minimizing the estimation error in the parameters of 
interest. Tuning parameter selection is performed using a 
multi-variable iterative search routine which seeks to 
minimize the theoretical mean-squared estimation error. This 
paper derives theoretical Kalman filter estimation error bias 
and variance values at steady-state operating conditions, and 
presents the tuner selection routine applied to minimize these 
values. Results from the application of the technique to an 
aircraft engine simulation are presented and compared to the 
conventional approach of tuner selection. Experimental 
simulation results are found to be in agreement with 
theoretical predictions. The new methodology is shown to 
yield a significant improvement in on-line engine performance 
estimation accuracy. 

Introduction 

An emerging approach in the field of aircraft engine 
controls and health management is the inclusion of real-time 
on-board models for the in-flight estimation of engine 
performance variations (Refs. 1, 2, and 3). This technology, 
typically based on Kalman filter concepts, enables the 
estimation of unmeasured engine performance parameters 
which can be directly utilized by controls, prognostics and 
health management applications. A challenge which 
complicates this practice is the fact that an aircraft engine’s 
performance is affected by its level of degradation, generally 
described in terms of unmeasurable health parameters such as 
efficiencies and flow capacities related to each major engine 
module. Through Kalman filter-based estimation techniques, 
the level of engine performance degradation can be estimated, 
given that there are at least as many sensors as parameters to 
be estimated. However, in an aircraft engine the number of 
sensors available is typically less than the number of health 
parameters presenting an under-determined estimation 
problem. A common approach to address this shortcoming is 


to estimate a sub-set of the health parameters, referred to as 
model tuning parameters. While this approach enables on-line 
Kalman filter-based estimation, it can result in “smearing” the 
effects of un-estimated health parameters onto those which are 
estimated, and in turn introduce error in the accuracy of 
overall model-based performance estimation applications. 
Recently, Litt (Ref. 4) presented an approach based on 
singular value decomposition which selects a model tuning 
parameter vector of low-enough dimension to be estimated by 
a Kalman filter. The model tuning parameter vector, q, was 
constructed as a linear combination of all health parameters, h, 
given by 

q = fh (1) 

where the transformation matrix, V , is selected applying 
singular value decomposition to capture the overall effect of 
the larger set of health parameters on the engine variables as 
closely as possible in the least squares sense. In this paper a 
new linear point design technique which applies a systematic 
approach to optimal tuning parameter selection will be 
presented. This technique, like the one presented in 
Reference 4, also defines a transformation matrix, V, used to 
construct a tuning parameter vector which is a linear 
combination of all health parameters, and of low enough 
dimension to enable Kalman filter estimation. The new 
approach optimally selects the transformation matrix, V , to 
minimize the theoretical steady-state estimation error in the 
engine performance parameters of interest. There is no known 
closed form solution for optimally selecting V to satisfy this 
objective. Therefore, a multivariable iterative search routine is 
applied to perform this function. 

The remaining sections of this paper are organized as 
follows. First, the mathematical formulation of the parameter 
estimation problem is presented, and theoretical estimation 
error values are derived assuming linear, steady-state 
operating conditions. The theoretical estimation error 
information is directly used by the iterative search routine 
applied to optimally select the Kalman filter tuning parameter 
vector, which is described next. Example estimation results 
from the application of the new methodology to an aircraft 
turbofan engine simulation are then presented and compared 
to the conventional approach of tuning parameter selection, 
the SVD tuner selection approach presented in Reference 4, 
and the maximum a posteriori performance estimation 
approach commonly applied for off-line (ground-based) 
aircraft engine gas path analysis applications (Refs. 5 and 6). 
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After the example, there is a discussion of practical 
considerations for applying the method, and a discussion of 
future work. Finally, conclusions are presented. 


Nomenclature 

A, A xh ,A xq 

A ■ B . /, . B , . 

6\/;, C x q, 

D, F, F xh , F xq , G, 
L,M,N 

system matrices 

C-MAPSS 

Commercial Modular Aero-Propulsion 
System Simulation 

Fn 

net thrust 

G x ■/;, G , . ( J , . Cl 

estimation bias matrices 

H 

matrix which relates health parameter 
effects to steady-state engine outputs 

HPC 

high pressure compressor 

HPT 

high pressure turbine 

/ 

identity matrix 

K x 

Kalman filter gain 

LPC 

low pressure compressor 

LPT 

low pressure turbine 

MAP 

maximum a posteriori 

Nf 

fan speed 

Nc 

core speed 

P h , r\ 

health & auxiliary parameter covariance 
matrices 

P24 

HPC inlet total pressure 

Ps30 

HPC exit static pressure 

^ih,k-’^,k,P:,k 

covariance matrices of estimated 
parameters 

P,. 

Kalman filter state estimation covariance 
matrix 

Q> Qxhi Qxq 

process noise covariance matrices 

R 

measurement noise covariance matrix 

SmLPC 

LPC stall margin 

T24 

HPC inlet total temperature 

T30 

HPC exit total temperature 

T40 

combustor exit temperature 

T48 

exhaust gas temperature 

T50 

LPT exit temperature 

V 

transformation matrix relating /?/ to q k 

vsv 

variable stator vane 

VBV 

variable bleed valve 

Wf 

fuel flow 

w z 

auxiliary parameter weighting matrix 

h k 

health parameter vector 



Kalman filter tuning parameter vector 

u k 

actuator command vector 

v k 

measurement noise vector 

Wk, w h ,k, W xh ,k 

process noise vectors 

Xk 

state vector 

Xxh,k 

augmented state vector (x k and /?;) 

%xq,k 

reduced order state vector (x t and (// ) 

yk 

vector of measured outputs 

Zk 

vector of unmeasured (auxiliary) outputs 

&xq,k 

residual vector (estimate minus its 

expected value) 

Subscripts 

k 

discrete time step index 

xh 

augmented state vector (x and h) 

xq 

reduced order state vector ( x and q) 

ss 

steady-state value 

Superscripts 

t 

pseudo-inverse 

A 

estimated value 

~ 

error value 

- 

mean value 

T 

transpose 

Operators 

E[-] 

expected value of argument 

H •} 

trace of matrix 

SSEE(-) 

sum of squared estimation errors 

WSSEE(-) 

weighted sum of squared estimation 

!! ■ IIf 

errors 

matrix Frobenius norm 


Problem Formulation 

The discrete linear time-invariant engine state space 
equations about a linear design point are given as 

x k +\ — Ax k + Bu k + Lh k + w k 

yk ~ C x k + D u k + Mh k + v k (2) 

z k = Fx k + Gu k + Nh k 

where k is the time index, x is the vector of state variables, u is 
the vector of control inputs, y is the vector of measured 
outputs, and z is the vector of auxiliary (unmeasured) model 
outputs. The vector h represents the engine health parameters, 
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which induce shifts in other variables as the health parameters 
deviate from their nominal values. The vectors w and v are 
uncorrelated zero-mean white noise input sequences. Q will be 
used to denote the covariance of w, and R to denote the 
covariance of v. The matrices A, B, C, D, F, G, L, M, and N 
are of appropriate dimension. The health parameters, 
represented by the vector h, are unknown inputs to the system. 
They may be treated as a set of biases, and are thus modeled 
without dynamics. With this interpretation Eq. (2) can be 
written as: 


x k + 1 


" A 

L~ 

x k 

+ 

~B~ 

u k + 

w k 

_hk + u 


0 

I 

k _ 


0 


_ Wh,k _ 



- A xh x xh,k + B xh u k + w xh,k 

+ D u k + v k 


y k = [c m] 

C xh 


x k 

h k 


~ C x h x xh,k + B>U k + V k 

+ Gu k 


z, = [F N] 


x k 

_ 

% xh ,k 

— ^ xh X xh ,k Gllk 


(3) 


The vector w xk is zero-mean white noise associated with the 
augmented state vector, [x T /? T ] T , with a covariance of Q xk . w xk 
consists of the original state process noise, w, concatenated 
with the process noise associated with the health parameter 
vector, Wh- 


^'xh , k 


W k 

Wh,k 


(4) 


The eigenvalues of A xk consist of the original eigenvalues of 
A plus an additional dim(/?) eigenvalues located at 1.0 on the 
unit circle due to the augmentation. Thus, the new augmented 
system given in Eq. (3) has at least as many eigenvalues 
located on the unit circle as there are elements of h. Once the h 
vector is appended to the state vector, it may be directly 
estimated, provided that the realization in Eq. (3) is 
observable. Using this formulation, the number of health 
parameters that can be estimated is limited to the number of 
sensors, the dimension of y (Ref. 7). Since in an aircraft gas 
turbine engine there are usually fewer sensors than health 
parameters, the problem becomes one of choosing the best set 
of tuners for the application. This paper presents a systematic 
methodology for the optimal selection of a model tuning 
parameter vector, q, of low-enough dimension to be estimated 
by a Kalman filter, while minimizing the estimation error in 
the model variables of interest. The following sub-sections 
will cover the steps in the problem setup. This includes 


construction of the reduced-order state space model, 
formulation of the Kalman filter estimator, calculation of the 
mean sum of squared estimation errors, and optimal selection 
of the transformation matrix to minimize the estimation error. 


Reduced-Order State Space Model 

The model tuning parameter vector, q, is constructed as a 
linear combination of all health parameters, h, given by 

q = V*h (5) 


where q e R m , h e IF, m < p, and V* is an m x p 
transformation matrix of rank m, applied to construct the 
tuning parameter vector. An approximation of the health 
parameter vector, h , can be obtained as 

h = V* T q (6) 

where V ' is the pseudo-inverse of V . Substituting Eq. (6) 
into Eq. (3) yields the following reduced order state space 
equations which will be used to formulate the Kalman filter 


x k + 1 


~A 

LV* t" 

x k 

3k+ 1_ 


0 

I 

3k _ 


u k + 


A xq Xxq,k 

— A X qX X q k + B X qU k + W X q k 

X k 


y k =[c MV* t' 


3k. 

*xq,k 


c m 

— C X qX X q k + Dll k + V k 


z k =[_F NV* t" 


x k 
3k _ 

X xq,k 


+ Gufc 


~ F X qXxq,k + GUfc 


w k 

Wq,k_ 

™ X q,k 


+ Dufc + 


(7) 


The state process noise, w xq , and its associated covariance, Q xq , 
for the reduced order system are calculated as 



7 

0 ' 


7 

o " 

w k 

II 

£ 

0 

V* 

W xh,k 

0 

V* 

3’h,k_ 


7 o " 


7 o " 


Qxh 


0 V* 

° v* _ 


Kalman Filter Formulation 

In this study, steady-state Kalman filtering is applied. This 
means that while the Kalman filter is a dynamic system, the 
state estimation error covariance matrix and the Kalman gain 
matrix are invariant — instead of updating these matrices each 
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time step they are held constant. Given the reduced order 
linear state space equations shown in Eq. (7), the state 
estimation error covariance matrix, P x , is calculated by 
solving the following Ricatti equation (Ref. 8): 


■Poo = AqPoo^xq - A xqPf.CL (C xq P x C^ q + R) 1 ... 


C X qPco A xq + Q 


(9) 


Cxq 


The steady-state Kalman filter gain, K., , can then be calculated 
as follows (Ref. 8) 

K x = P x C T xq {C xq P y fi T xq + R)' (10) 

and, assuming steady-state, open-loop operation (u = 0), the 
Kalman filter estimator takes the following form 

X xq,k ~ A xq x xq,k—l + ( }'k ~ C X q A xq x xq,k-l ) (H) 


The reduced order state vector estimate, x xq , produced by 

Eq. (11) can be used to produce an estimate of the augmented 
state vector, and the auxiliary parameter vector as follows 


X xh,k 


I 

0 



z k=[ E NV*^x^ k 


( 12 ) 


Analytical Derivation of Estimation Error 

The estimation errors in x x h£ and zy. are defined as the 
difference between the estimated and actual values) 


where the operator £[•] represents the expected value of the 
argument. The variance of the estimates can be found by 
constructing their respective estimation covariance matrices 


P ih.k t 


P~z,k=E\ 


( X xh,k - E \xxh,k]){xxh,k - E [xxh,k ]) 7 

(zk- E [h]){h- E [zk ]f 


(15) 


Diagonal elements of the covariance matrices will reflect the 
variance in individual parameter estimates, while off-diagonal 
elements reflect the covariance between parameter estimates. 
The overall sum of squared estimation errors (SSEE) can be 
obtained by combining the estimation error bias and 
estimation variance information as 

SSEE (i ' x h,k ) = X x h,k*xh,k + tr \ E kh,k } ^ ^ 

SSEE(z k ) = z[z k +tr{P= Jc } 

where tr { . j represents the trace (sum of the diagonal 
elements) of the matrix. In this paper, theoretical values for 
each error component will be derived assuming steady-state, 
open- loop (u = 0) operating conditions. First, the estimation 
error bias is derived, followed by a derivation of the 
estimation variance. 

Estimation Error Bias 

The estimation error biases, 'x x h,k an d z k , can be 

analytically derived for an arbitrary health parameter vector, h, 
at steady-state operating conditions. This is done taking 
advantage of the following expected value properties at 
steady-state open-loop operating conditions 


X xh,k ~ X xh,k X xh,k 
% = Z k ~ z k 


(13) 


Due to the under-determined nature of the estimation problem, 
it will be impossible for the Kalman filter estimator to 
completely restore all information when transforming q into 

h . As such, the Kalman filter will be a biased estimator (i.e. 
the expected values of x x /, ^ an ^ z"/. will be non-zero). The 

estimation errors can be considered to consist of two 
components: An estimation error bias, and an estimation 
variance. The estimation error bias vectors are equivalent to 
the mean estimation error vectors defined as 

X xh,k ~ E \ X xh.k \ — E \ X xh,k ~ X xh,k 

, (14) 

z k = E [ z k] = E [ z k - z k ] 


E [-'■i +l ] — E [ x k ] — x ss 

E[h] = h 

E \_ x xh,k ] — X xh,ss 

E [}’k] = y ss 

E [ z k] = z ss 
E [u k ] = 0 
E[w k ] = 0 
E [v k ] = 0 


E \_ x xq.k J E \_ x xq,k-\ J x . 
E \ x xh,k ] — X xh,ss 


xq,ss 


E [ z k] = z 


ss 


(17) 
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where the subscript “ss” denotes steady-state operation. By 
taking expected values of Eq. (2), x ss , y ss and z ss can be written 
as functions of the health parameter vector h 


E[*k+i\ = A- E[x k ] + B ■ E[u k ] + L ■ E[h k ] + E[w k ] 

x ss = Ax ss + Lh 

x ss ={l-Ay l Lh 


(18) 


E[yk] = C-E[x k ] + D-E[u k ] + M-E[h k \ + E[v k ] 

yss ~ Cx s s Mh 

y ss =(c{l-A)~ l L+My 


E[z k ] = F ■ E [x k ] + G ■ E \u k \ + N ■ E[h k ] 
z ss = Fx ss + Nh 

z ss =[f(i-a)- 1 l+nY 


Next, by taking expected values of both sides of Eq. (11), the 
expected value of x xq ^ can be obtained as a function ofy ss : 


F X xq,k ^xq ' F^X xq , k —\ J • • • 


*xq,ss 
% xq,ss 


+ ^oo ] C X qA X q ’ j 

— A X qX X q SS + K„ ^55 — C X q-AxqX X q,ss j 

— — A X q + C X q A X q ^ K^y ss 


( 21 ) 


Then, making the substitution = (c(/ -A)~ X L + Mji given 
in Eq. (19), the expected steady-state value of x xq ^ can be 
written as a function of h 


% xq,ss 


A X q + C X q A X q ^ 

X kAc(I - A)~ l L + M\h 


(22) 


The steady-state augmented state estimation error bias can 
then be found, and partitioned into error bias information for 

the original state vector, x ss , and the health parameter vector, 
h ss , by combining Eqs. (12), (14), (18) and (22) to yield 


X xh,. 

X xh,. 


F X xh,k X xh,k 


X xh,ss 


X xh,ss X xh,ss 
x ss 


X xh,ss 


n ss 

I 0 

0 F*t 


X xq,ss X xh,ss 


I 0 

0 F*t 


(23) 


P ^xq KazCxqAxq ) 


xK„ 


C(I-A) 1 L+M 


(I-A)~ { L 
I 

\r~ 

G x i, 


X xh,ss 


O x h h 


The steady-state auxiliary parameter estimation error bias can 
also be derived by combining Eqs. (12), (14), (20), and (22) to 
yield 


^SS — E\zfc z k 


^SS ^SS ^ss 


= [f NV* 

[f NV* t ] (/ - A xq + K x C xq A xq 
C[l-A]~' L+M 


xK a 


(24) 


F(I-A) 1 L+N 


z ss = G : h 


The estimation error bias equations given in Eqs. (23) and 
(24) are functions of an arbitrary health parameter vector, h. 
As such they are representative of the parameter estimation 
error biases in a single engine, at a given point in its lifetime 
of use where its deterioration is represented by the health 
parameter vector h. The average sum of squared estimation 
error biases across a fleet of engines can be calculated as 
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xh, fleet 


= E 


_*xh,ss*xh,ss _ E tr {xxh,ss*xh,ss} 


= E[tr{G xh hhTGl h }] 
G xh ■ E[hh T ]G T xh 


= tr< 


(25) 


= tr{G xh P h G T xh 


6 fleet 


— E jz^z^ } J 

= E[tr{G z hh T G T z ]\ 


= tr{ 


G z -E[hh T ]Gl 

"~p h 7 ' 

= tr{G.P h Gl } 


(26) 


where the matrix P h , defined as E[hh‘ ], reflects a priori or 
historical knowledge of the covariance in the health 
parameters across all engines. If available, it can be used to 
predict the sum of squared estimation errors biases as shown 
in (25) and (26). 

Estimation Variance 

Next, derivations are presented for the augmented state 
estimate and auxiliary parameter estimate covariance matrices, 
P- h k and P~ zk respectively. These matrices will be 

calculated as a function of the reduced-order state vector 
estimation covariance matrix, P xq ,k > which is defined as 


E’xq,k ~ E 


(x X q,k -E[x xq ^[x xq ^ k -E\jc xq ,k^) 


E-xq.k 


f T 

°xq,k 


(27) 


where the vector s xqk is defined as the residual between x xqk 

at time k and its expected value. Since E\x xq k ] = x xqss , z xqk 
can be obtained by subtracting Eq. (21) from Eq. (11) 


E\jX X q,k \ 

- ^xq,SS 

— A X qX X q k- 1 + ^oo ( yk ~ C X q-A X qX X q,k- 1 ) 


£ X q,k X X qk 

~ % xq,k — -^xq,ss 


(28) 


*xq,k 

~ {^xq^xq,ss + ^oo ~ C X qA X qX X q SS 

K 3~V 

%xq,ss 

~ ( ^xq ~ K x CxqA X q X xqk _\ — X xqxx j 

+K< x (y k -y ss ) 


Making the substitutions z xq ,k-i = *xq,k-i ~*xq,ss > and 

n =yk -y ss yields 


'-'xq,! 


(30) 


i,k iy^xq KtoCxqAxq^Zxq^k—X+K-coVk (29) 
The estimation covariance matrix P X q,k is then calculated as 

E X q,k = 

= \^xq ~ KaoCxqAxq 

X [ ^xq ~ K x C X q A X q J ... 

\_^xq ~ EaoCxqA X q J E [Twji.i’-A/t J E x . . 

+ K x E\v k z T xq ^[A xq K(x> C X q A X q J ... 

+ K„E[v k v T k ]Kl 

The substitutions E^ xq ,k-^ T xq , k -i\= p iq,k-\, and E\y k vl\=R 
can be made in the above equation. Since £ xq , k -\ and are 

uncorrelated, the substitution 7s[ s *?,*-i v ifc ] = d v * s x?,/t-i] = ® 
can also be made, producing 

lT 


E x q,k A X q K x C xq A xq P xq ^k-\ ^xq E Xj C xq ^xq J 


+ K x RKl 


(31) 


At steady-state operating conditions P xq ,k- 1 = Exq,k ■ Making 
this substitution in (31) produces the following Ricatti 
equation which can be solved for P kqk ■ 


E\q , k ^xq E x C xq A xq P xq k A xq K 0 0 C xq A- 


xq J 


+ k,rkI 


(32) 
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It should be noted that P xqk obtained by solving (32) will be 

identical to P x produced via Eq. (9) if the system’s actual state 
process noise covariance is identical to the Q xq assumed in the 
design of the Kalman filter. However, Q is often treated as a 
Kalman filter design parameter to provide acceptable dynamic 
response. For the purpose of this derivation, we have assumed 
a steady-state operating condition where the state variables 
and health parameters are invariant, and thus the actual system 
process noise is zero (i.e., w xk ,k= 0)- In this case P xqk will not 

equal P x . Once P X q^k * s obtained, it can be used to calculate 
P~l , , the covariance of x xkk , which is defined as 



7 

o " 

p , 

7 o " 

0 

v*t_ 

xq,k 

i 

O 

l 


(33) 


The augmented state vector estimation covariance given in Eq. 
(33) can be partitioned into covariance information for the 
original state vector, P% k (upper left corner of the P~; , 

matrix), and the health parameter vector, P; (lower right 

corner of the P-; , matrix) 


SSEE {x xk fl ee t ) + tr \^xh,k } 

= tr { G xhPhG T xh +P iilk } (36) 
SSEE(z aeet ) = ^l eet +tr{P^ k } 

= tr{G z P h Gj+P^ 4 

If required, a weighted sum approach can be applied to 
normalize the contributions of individual auxiliary parameter 
estimation errors. This is often necessary as there may be 
several orders of magnitude difference between the auxiliary 
parameters of interest. A weighted sum approach prevents 
domination by individual parameters. In this study a diagonal 
auxiliary parameter weighting matrix, W z , is applied based on 
the inverse of auxiliary parameter variance (obtained from the 
main diagonal of the auxiliary parameter covariance matrix, 

P : ) 


P = 


W- = 


F(I-A) 1 L + N P h F(I-A) 1 L + N 


~\T 


Pz, it 0 0 

0 0 
0 0 P 7h 


-1 -1 


(37) 


^xh.k 



(34) 


W : is then applied to calculate of a “weighted” sum of 
auxiliary parameter squared estimation errors given as 


The P’ X q.k matrix from Eq. (32) can also be used to calculate 
Pz k , the covariance in the estimation of z k , which is 

equivalent to -E[z k ]\zk 

Pz,k=[F NV*^P^ k [F NV* t] r (35) 

The variance in the estimates x xh k and z k can be obtained 

from the diagonals of the covariance matrices produced by 
(33) and (35) respectively. 

Sum of Squared Estimation Errors 

Once Eqs. (25), (26), (33) and (35) are obtained, they may 
be used to analytically calculate the mean sum of squared 
estimation errors over all engines by combining the respective 
estimation error bias and estimation variance information as 
previously shown in Eq. (16). The mean augmented state 
vector sum of squared estimation errors, SSEe(x x/i (lee[ ) , and 
the mean auxiliary parameter vector sum of squared 
estimation errors, .S',S7:f (z ncct ), become 


WSSEE (zfl eet ) = tr{w z [ G : P h GT + P k ]} (38) 

From Eqs. (23), (24), (33) and (35) it can be observed that 
both bias and variance are affected by the selection of the 
transformation matrix, V . The sum of squared estimation 
error terms derived in this section give rise to an optimization 
problem: selecting V to minimize the squared estimation error 
in the Kalman filter produced parameter estimates. This could 
include health parameter estimates, auxiliary parameter 
estimates, or a combination of parameters. Although there is 
no known closed form solution for optimally selecting the V * 
matrix to satisfy the objective of minimizing estimation errors, 
a multi-parameter iterative search method has been developed 
to perform this task, and will be described in the next section. 

Optimal Transformation Matrix Selection 

Prior to initiating the search for an optimal V* , specific 
system design information must be defined or obtained. This 
includes: 

• Specifying the auxiliary parameters to be estimated 

• Generating system state space equations at a fleet average 
(50 percent deteriorated) engine trim point 
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• Defining measurement noise covariance matrix, R 

• Defining augmented state process noise covariance 
matrix, Q xh 

• Defining fleet average health parameter covariance, P h 

Some additional clarification is provided regarding the 
selection of Pi, and Q x h as the distinction between these two 
covariance matrices may not be immediately obvious. Pi, 
defines the expected health parameter covariance across all 
engines. It may be based on past knowledge gained from 
engine gas path analysis programs and/or historical studies of 
engine module performance deterioration. Conversely, Q x t, 
defines the expected process noise covariance in the state 
variables and health parameters of an individual engine, at a 
single discrete time step, k. The selection of Q xh will directly 
impact the dynamic response and the variance of the estimates 
generated by the Kalman filter, and to a large extent Q x /, is 
treated as a design parameter. 

After the necessary system information has been obtained, 
the search for an optimal transformation matrix to minimize 
the Kalman filter sum of squared estimation errors can 
commence. This is performed using the Isqnonlin function of 
the Matlab(The MathWorks, Inc.) Optimization Toolbox. This 
function applies an iterative search to find the least squares 
solution of a user-specified multivariable optimization 
problem. A flow chart depicting the steps in this optimal 
iterative search is shown in Figure 1, and a further description 
of each step is given below. 

1. Upon startup, an initial random guess of V is 
generated. It is selected such that the matrix Frobenius 

norm \ F = 1. This requirement is applied to help 

prevent convergence to a poorly scaled solution. 

2. Construct reduced order state-space model, (Eq. (7)). 

3. Formulate Kalman filter 

a. Calculate estimation covariance matrix, P w (Eq. (9)) 

b. Calculate the Kalman gain matrix, K*,, (Eq. (10)) 

4. Calculate sum of squared estimation errors (Eq. (36)), or 
weighted sum of squared estimation errors (Eq. (38)) 

5. On each iteration the change in SSEE (or WSSEE) 
relative to the previous iteration is assessed to 
determine if convergence within a user specified 
tolerance has been achieved. 

a. If converged, skip step 6 and proceed to step 7 

b. If not converged, proceed to step 6 to update V 

6. V is updated via the Matlab Isqnonlin function, again 

requiring that ||^ || = 1, an d the process returns to 
step 2 

7. Upon convergence, the optimization routine returns 
the optimal value of V , and ends. 

Experience has shown that the transformation matrix 
returned by the optimization routine is not unique — different 
matrices can be found which produce a global minimum of the 



Figure 1 . — Flowchart of V Iterative Optimal Search 


objective function. Experience has also shown that the 
optimization routine will usually return a V matrix which 
satisfies, or nearly satisfies (i.e., within 5 percent), the global 
minimum of the objective function. However, in order to 
guard against potential convergence to a local minimum, it is 
prudent to run the optimization routine multiple times, each 
time starting with a different initial guess for V . This is only 
to assure the designer that the global minimum is achieved, 
not to produce a consistent Y . It should be emphasized that 
the optimal search for Y is only conducted off-line during the 
estimator design process. This calculation is not conducted as 
part of the on-line real-time Kalman filter implementation, and 
thus places no additional computational burden upon it. 

Turbofan Engine Example 

A linearized cruise operating point extracted from the 
NASA Commercial Modular Aero-Propulsion System 
Simulation (C-MAPSS) high-bypass turbofan engine model 
(Ref. 9) is used to evaluate the new systematic tuner selection 
methodology. The linear model has two state variables, ten 
health parameters, and three control inputs, all shown in 
Table 1. The model’s seven sensed outputs, and corresponding 
sensor noise standard deviation, are shown in Table 2. The 
auxiliary output parameters of interest to be estimated are 
shown in 

Table 3. The linear model is used as the truth model for this 
application example. The model is run open-loop, so all 
control inputs remain at 0, i.e., they do not deviate from the 
trim value for the linear model and no actuator bias is present. 
Deviations in all ten health parameters are assumed to be 
uncorrelated, and randomly shifted from their trim conditions 
with a standard deviation of ±0.02 (±2 percent). Since a 
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parameter’s variance is equal to its standard deviation squared, 
the health parameter covariance matrix, P/„ is defined as a 
diagonal matrix with all diagonal elements equal to 0.0004. 

Next, the estimation accuracy of the systematic approach 
for selecting Kalman filter tuning parameters will be 
compared to the conventional approach of selecting a sub-set 
of health parameters to serve as tuners (the seven health 
parameters denoted with in Table 1), and the singular 
value decomposition approach to tuner selection introduced by 
Litt in Reference 4. Table 4 shows a comparison of the 
theoretically predicted estimation errors (squared bias, 
variance, and total squared error) and experimentally obtained 
squared estimation errors for each of the three tuner selection 
approaches. T40 and T50 estimation errors are shown in 
squared degrees Rankine, and Fn and SmLPC estimation 
errors are shown in squared percent net thrust and squared 
percent stall margin respectively. The experimental results 
were obtained through a Monte Carlo simulation analysis 
where the health parameters varied over a random distribution 
in accordance with the covariance matrix, P/,. The test cases 
were concatenated to produce a single time history input 
which was provided to the C-MAPSS linear discrete state 
space model given in Eq. (2), with an update rate of 15 ms. 
Each individual health parameter test case lasted 30 s. 


TABLE 1.— STATE VARIABLES, HEALTH 
PARAMETERS, AND ACTUATORS 


State variables 

Health parameters 

Actuators 

Nf — fan speed 
Nc — core speed 

Fan efficiency 
Fan flow capacity* 
LPC efficiency* 
LPC flow capacity 
HPC efficiency* 
HPC flow capacity* 
HPT efficiency* 
HPT flow capacity* 
LPT efficiency 
LPT flow capacity* 

Wf — fuel flow 

VSV — variable stator vane 

VBV — variable bleed valve 


*Health parameters selected as tuners in conventional estimation approach, 
(identical to sub-set of health parameters selected as tuners in Ref. 4) 


TABLE 2.— SENSED OUTPUTS AND STANDARD DEVIATION 
AS PERCENT OF OPERATING POINT TRIM VALUES 


Sensed output Standard deviation 

(%i 

Nf- — fan speed 0.25 

Nc — core speed 0.25 

P24 — HPC inlet total pressure 0.50 

T24 — HPC inlet total temperature 0.75 

Ps30 — HPC exit static pressure 0.50 

T30 — HPC exit total temperature 0.75 

T48 — Exhaust gas temperature 0.75 


TABLE 3.— ESTIMATED AUXILIARY 
PARAMETERS 

Auxiliary parameter 

T40 — Combustor exit temperature 
T50 — LPT exit temperature 
Fn — Net thrust 
SmLPC — LPC stall margin 


At the completion of each 30 s test case, the health 
parameter vector input instantaneously transitioned to the next 
test case. A total of 375 30 s test cases were evaluated, 
resulting in an 11,250 s input time history. Three separate 
Kalman fdters were implemented using the three tuner 
selection approaches. The experimental estimation errors were 
determined by calculating the mean squared error between 
estimated and actual values during the last 10 s of each 30 s 
test case. The error calculation is based on only the last 10 s so 
that engine model outputs and Kalman estimator outputs have 
reached a quasi-steady-state operating condition prior to 
calculating the error. This ensures that the experimental results 
are consistent with the theoretically predicted estimation errors 
which were derived assuming steady-state operation. 

From Table 4 it can be seen that the theoretically predicted 
and the experimentally obtained squared estimation errors 
exhibit good agreement. If the number of random test cases 
were increased to a suitably large number, it is expected that 
the theoretical and experimental results would be identical. It 
can also be seen that all three estimators are able to produce 
unbiased estimates of the combustor exit temperature, T40; 
however, their estimates of LPT exit temperature, T50, net 
thrust, Fn, and LPC stall margin, SmLPC, are biased. The 
encouraging finding is that the new systematic approach to 
tuner selection significantly reduces the overall mean squared 
estimation error compared to the other two approaches. 
Relative to the conventional approach of tuner selection the 
experimental mean squared estimation errors in T40, T50, Fn 
and SmLPC are reduced 76, 82, 80 and 63 percent, 
respectively. It can also be observed that the SVD tuner 
selection approach, which is designed to reduce the estimation 
error bias, does in fact reduce the sum of squared biases 
relative to the subset of health parameters approach. However, 
the SVD approach is also found to increase the estimation 
variance, which contributes to its overall mean squared 
estimation error. 


TABLE 4.— AUXILIARY PARAMETER 
SQUARED ESTIMATION ERRORS 


Tuners 

Error 

T40 

T50 

Fn 

SmLPC 



(°R) 

(°R) 

(%) 

(%) 

Subset of 
health 
parameters 

Theor. sqr. bias 

0.00 

561.76 

3.84 

3.28 

Theor. variance 

74.76 

29.65 

0.48 

0.34 

Theor. sqr. error 

74.76 

591.41 

4.31 

3.62 

Exper. sqr. error 

74.90 

583.29 

4.27 

3.60 


Theor. sqr. bias 

0.00 

512.46 

4.05 

5.28 

SVD tuner 

Theor. variance 

65.99 

67.21 

0.80 

1.31 

selection 

Theor. sqr. error 

65.99 

579.67 

4.86 

6.59 


Exper. sqr. error 

66.20 

579.39 

4.98 

6.76 

Systematic 

tuner 

selection 

Theor. sqr. bias 

0.00 

87.81 

0.66 

0.95 

Theor. variance 

17.49 

18.55 

0.13 

0.35 

Theor. sqr. error 

17.49 

106.35 

0.79 

1.30 

Exper. sqr. error 

17.61 

106.54 

0.86 

1.35 


A visual illustration of the effect that tuner selection has on 
Kalman filter estimation accuracy can be seen in Figures 2 
through 5, which show actual and estimated results for the 
auxiliary parameters T40, T50, Fn and SmLPC respectively. 
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Each plot shows a 300 s segment of the evaluated test cases. 
The step changes that can be observed in each plot every 30 s 
correspond to a transition to a different health parameter 
vector. True model auxiliary parameter outputs are shown in 
black, and Kalman fdter estimates are shown in red. In each 
figure the information is arranged top to bottom according to 
tuner selection based upon: a) a subset of health parameters; b) 
singular value decomposition; and c) the new systematic 
selection strategy. The information shown in these figures 
corroborates the infonnation in Table 4; namely all three tuner 
selection approaches produce unbiased estimates of T40 
(Figure 2), while the systematic tuner selection strategy yields 
a noticeable reduction in the total squared estimation error 
(squared bias plus variance) of all four auxiliary parameters. 

Comparison with Maximum a Posteriori 
Estimation 

The presented systematic tuner selection strategy minimizes 
the mean squared error of the on-line estimator at steady-state 
operating conditions, taking advantage of prior knowledge of 
engine health parameter distributions. As such it is somewhat 
analogous to the maximum a posteriori (MAP) estimation 
method which is commonly applied for ground-based aircraft 
gas turbine engine gas path analysis (Refs. 5 and 6). This leads 
to the question, how does the on-line Kalman filter estimation 
accuracy compare to MAP estimation accuracy? Prior to 
making this comparison the mathematical formulation of the 
MAP estimator is briefly introduced. Here a steady-state 
model of the measurement process in the following form is 
applied 


results. To test this theory, a MAP estimator was designed and 
its estimation accuracy was compared to a Kalman filter with 
tuning parameters optimally selected to minimize the 
estimation errors in the health parameter vector h. First, the 
two estimators were designed and evaluated using the original 
sensor noise levels shown in Table 2. Next, the sensor noise 
levels were set to l/20 ,h of their original levels, the estimators 
were re-designed, and the comparison was repeated. Monte 
Carlo simulation evaluations as previously described were 
applied (i.e., 375 random health parameter vectors, 30 s in 
duration, with estimation accuracy calculations based upon the 
last 10 s of each 30 s test case). Theoretical and experimental 
estimation errors are shown in Table 5 and Table 6 for the 
original noise and reduced noise levels, respectively. At 
original noise levels the Kalman estimator is able to produce 
smaller estimation errors. However, at the reduced noise level 
the two estimation approaches are found to be nearly identical. 
This comparison validates that the Kalman estimation 
approach is indeed producing a minimum mean squared 
estimation error as intended, while providing the capability to 
support real-time on-line estimation under dynamic operating 
scenarios. 


Subset of health parameters 


cr 50 
< -50 
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- Actual 
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SVD tuner selection 
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< -50 
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yk - Hh k + v k (39) 

where the matrix H relates the effects of the health parameter 
vector, h, to the sensed measurements, y. From Eq. (19), it can 
be seen that H is equivalent to C(1 — A)~ 1 L + M. The maximum 
a posteriori (MAP) estimator follows the closed form 
expression 


h k =[p^+H t R~ x h) 1 H T R~ l y k (40) 

The MAP estimator is capable of estimating more 
unknowns than available measurements due to the inclusion of 
a priori knowledge of the estimated parameter covariance, P/,. 
However, the MAP estimator, unlike a Kalman filter, is not a 
recursive estimator and does not take advantage of past 
measurements to enhance its estimate at the current time step. 
Furthermore, the MAP estimator only considers a static 
relationship between system state variables and measured 
outputs — it does not consider system dynamics. Because of 
these differences a Kalman estimator with optimally selected 
tuning parameters should outperform the MAP estimator. 
However, under steady-state conditions, with minimal sensor 
noise the two estimation approaches should produce similar 


Systematic tuner selection 
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Figure 2. — T40 estimation (tuner comparison). 


Subset of health parameters 


- Estimate 
Actual 


A | - 

-50 - Lyyii 

0 50 100 150 200 250 300 


SVD tuner selection 


0 50 100 150 200 250 300 


Systematic tuner selection 



0 50 100 150 200 250 300 


time (sec) 

Figure 3. — T50 estimation (tuner comparison). 
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Figure 4. — Fn estimation (tuner comparison). 


Figure 5. — SmLPC estimation (tuner comparison). 


TABLE 5.— HEALTH PARAMETER PERCENT SQUARED ESTIMATION ERRORS (NOMINAL NOISE) 


Estimator 

Error type 

hi 

h2 

hi 

h4 

h5 

h6 

h7 

h8 

h9 

hlO 

Sum 


Theor. squared bias 

1.26 

1.48 

2.28 

1.22 

0.00 

0.00 

0.74 

0.00 

1.97 

3.06 

12.00 

Kalman 

Theor. variance 

0.73 

0.19 

0.80 

0.45 

0.40 

0.68 

0.21 

0.21 

0.34 

0.08 

4.09 

filter 

Theor. squared error 

2.00 

1.67 

3.07 

1.66 

0.40 

0.68 

0.96 

0.21 

2.32 

3.14 

16.09 


Exper. squared error 

1.82 

1.60 

2.77 

1.59 

0.41 

0.68 

0.97 

0.21 

2.13 

3.15 

15.31 


Theor. squared bias 

2.42 

1.75 

3.53 

1.90 

0.47 

0.96 

1.00 

0.16 

2.47 

3.17 

17.81 

MAP 

Theor. variance 

0.22 

0.27 

0.16 

0.64 

0.69 

0.71 

0.36 

0.36 

0.17 

0.09 

3.68 

estimator 

Theor. squared error 

2.63 

2.02 

3.69 

2.54 

1.16 

1.67 

1.36 

0.52 

2.64 

3.26 

21.48 


Exper. squared error 

2.52 

1.89 

3.44 

2.45 

1.14 

1.57 

1.34 

0.51 

2.49 

3.24 

20.60 


TABLE 6.— HEALTH PARAMETER PERCENT SQUARED ESTIMATION ERRORS (REDUCED NOISE) 


Estimator 

Error type 

hi 

h2 

h3 

h4 

h5 

h6 

h7 

h8 

h9 

hlO 

Sum 


Theor. squared bias 

1.26 

1.48 

2.28 

1.22 

0.00 

0.00 

0.74 

0.00 

1.97 

3.06 

12.00 

Kalman 

Theor. variance 

0.02 

0.00 

0.03 

0.01 

0.01 

0.02 

0.00 

0.00 

0.01 

0.00 

0.10 

filter 

Theor. squared error 

1.29 

1.48 

2.30 

1.22 

0.01 

0.02 

0.74 

0.00 

1.98 

3.06 

12.10 


Exper. squared error 

1.13 

1.41 

1.97 

1.16 

0.01 

0.02 

0.75 

0.00 

1.79 

3.08 

11.31 


Theor. squared bias 

1.26 

1.48 

2.28 

1.22 

0.00 

0.00 

0.74 

0.00 

1.97 

3.06 

12.00 

MAP 

Theor. variance 

0.03 

0.00 

0.04 

0.01 

0.01 

0.02 

0.00 

0.00 

0.01 

0.00 

0.13 

estimator 

Theor. squared error 

1.30 

1.48 

2.32 

1.22 

0.01 

0.02 

0.74 

0.00 

1.99 

3.06 

12.14 


Exper. squared error 

1.14 

1.41 

1.98 

1.17 

0.01 

0.02 

0.75 

0.00 

1.79 

3.08 

11.34 


Discussion 

While the systematic tuner selection approach presented 
here appears promising for on-line Kalman filter based 
parameter estimation applications, there are several practical 
considerations which need to be assessed when applying such 
a technique. The optimization routine attempts to minimize the 
overall squared estimation error — both bias and variance — 
under steady-state operating conditions. The minimization of 
the estimation variance in particular can come at the expense 
of dynamic responsiveness of the Kalman filter. To illustrate 
this consider the time history plots of actual versus estimated 
T40 shown in Figure 6. The top plot shows Kalman 
filter estimation results using a tuning parameter vector 


systematically selected to minimize the error in four auxiliary 
parameters (T40, T50, Fn and SmLPC) as presented in the 
previous section. The bottom plot shows Kalman filter 
estimation results using a tuning parameter vector 
systematically selected to minimize the estimation error in 
T40 only. At time 100 s a step change in the health parameter 
input vector is introduced into the engine model; this allows 
the dynamic response of the two estimators to be compared. 
It can be observed that T40 estimation variance in the 
bottom plot is reduced, as is the mean steady-state estimation 
error (> 300 s). This is not surprising since one would 
generally expect improved results when optimizing 
to minimize the error in a single parameter, as opposed 
to multiple parameters. However, the estimator shown in the 
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Figure 6. — Illustration of tuner impact 
on estimator response. 

bottom plot does require a significantly longer time to reach 
steady-state convergence. Conversely, the estimator designed 
to minimize the steady-state error in four auxiliary parameters 
(top plot) is unable to place as much emphasis on T40 
estimation variance reduction, but it is able to track dynamic 
changes in T40 more rapidly. This example illustrates the 
inter-dependence between estimation variance and 
responsiveness. Therefore, it is prudent for a designer to 
evaluate the Kalman filter to ensure that it tracks engine 
dynamics acceptably. If the dynamic response is unacceptable, 
the optimization routine can be rerun placing more weight on 
estimation error bias reduction, and less weight on variance 
reduction. 

All results presented in this paper are based on a turbofan 
engine linear state-space model at a single operating point. 
While this linear assessment enables experimental validation 
of theoretical predictions, it is not representative of actual 
aircraft gas turbine engine operation which is a transient non- 
linear system operating over a broad range of operating 
conditions. In order for the presented approach to be 
applicable it would need to be able to produce an optimal set 
of tuning parameters, not just at a single operating point but 
rather a globally optimal tuning parameter vector universally 
applicable over the range of operating conditions that an 
engine is expected to experience. A potential approach to 
selecting a single “globally optimal” tuning parameter vector 
is to modify the optimization routine to minimize the 
combined estimation error over multiple engine operating 
points such as takeoff, climb and cruise. This would be a 
straightforward modification to the Matlab optimization 
routine, but it would increase the computational time required 
to calculate the result. Since the systematic tuner selection 
process is only envisioned to be done once during the system 
design process, this will not impact the on-line execution 


speed of the Kalman filter. It is anticipated that the application 
of globally optimal tuners will result in some estimation 
accuracy degradation relative to tuners optimized for 
individual operating points, although this has not yet been 
verified or quantified. 

Conclusions 

A systematic approach to tuning parameter selection for on- 
line Kalman filter based parameter estimation has been 
presented. This technique is specifically applicable for the 
underdetermined aircraft engine parameter estimation case 
where there are fewer sensor measurements than unknown 
health parameters which will impact engine outputs. It creates 
and applies a linear transformation matrix, V, to select a 
vector of tuning parameters which are a linear combination of 
all health parameters. The tuning parameter vector is selected 
to be of low-enough dimension to be estimated, while 
minimizing the mean-squared error of Kalman filter estimates. 
The multi-parameter iterative search routine applied to 
optimally select V was presented. Results have shown that 
while the transformation matrix returned by the optimization 
routine is not unique (different matrices can be found which 
produce a global minimum of the objective function), the 
routine is effective in returning a transformation matrix which 
is optimal, or near optimal, regardless of its initial starting 
guess of the matrix. The efficacy of the systematic approach to 
tuning parameter selection was demonstrated by applying it to 
parameter estimation in an aircraft turbofan engine linear point 
model. It was found to significantly reduce mean squared 
estimation errors compared to the conventional approach of 
selecting a subset of health parameters to serve as tuners. In 
some parameters the mean squared estimation error reduction 
was found to be over 80 percent. These estimation 
improvements were theoretically predicted and experimentally 
validated through Monte Carlo simulation studies. 

The systematic approach to Kalman filter design is 
envisioned to be applicable for a broad range of on-board 
aircraft engine model-based applications which produce 
estimates of unmeasured parameters. This includes model- 
based controls, model-based diagnostics, and on-board life 
usage algorithms. It is also envisioned to have benefits for 
sensor selection during the engine design process, specifically 
for assessing the performance estimation accuracy benefits of 
different candidate sensor suites. Areas for future work 
include extending the technique to produce a tuning parameter 
vector optimal over a range of operating conditions, and 
evaluating the technique on a non-linear engine model, under 
both steady-state and transient operating conditions. 
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